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1. Background and Objectives 


• As the normal angle of attack, normal Mach number and sweep-back angle are 
varied; complex flows develop around Delta and belta-like wings. 

These flows are characterized with the formation of large and small scale 
vortices , weak and strong shock waves , and shock induced separations . The 
experimental Miller and Wood Classification Diagram shows seven regions of this 
flow. 


• These flows become highly complex when vortex breakdown occurs in the vicinity of 
the wing or when the wing undergoes unsteady motion due to maneuvering or 
flutter. 


• The main objectives of this ongoing research research work are to develop 

efficient and reliable computational schemes which are capable of predicting the 
distributed aerodynamic characteristics of these wings in steady and unsteady 
flows over a wide range of angles of attack, sweep-back angles, Mach numbers and 
configurations. 





1 Classical 
Vortex 

2 Separation 
Bubble with 
Mo Shock 

3 No Shock/ 

No Separation 

4 Shock with no 
Separation 

5 Shock- Induced 
Separation 

6 Separation 
Bubble with 
Shock 

7 Vortex with 
Shock 


Fig. l Miller and Wood 1 Classification Diagram. 
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2. Unsteady Euler Equations in a Rotating Frame of 
Reference for Transonic-Vortex Flows 

Formulation 

• Conservation Form of Euler Equations in a Space-Fixed Frame of Reference 


*i + &!+ o 

(1) 

at ax ay az 

where 

q = [p, pu, pv, pw, pe] 1 

(2) 

2 t 

F = [pu, pu + p, puv, puw, puh] 

(3) 

F = [pv, puv, p v 2 + p, pvw, pvh] t 

(4) 

2 _,t 

<a = [pw, puw, pvw, pw + p, pwhj 

(5) 

e = p/p(y-l) + (u 2 + v 2 + w 2 )/? 

(6) 

h = e + p/p 

(7) 

• Rewriting the Equations in the Vector Form 
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(9) 

M + V • (p h V) - 0 

(10) 
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And Using the Substantial and Local Derivatives Relations 


Da D'a aa a' a 
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( 12 ) 


where V 
Frame 


- - D' a' 

V + wxr ; — - , — = Substantial and Local Derivative in the Rotating 
r nt 


• We get the Conservative Form of Euler Equations for the Relative Motion 
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Method of Solution 


• The Abstract Conservative Form of the Relative Motion in Terms of Rotating 
Coordinates are 
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The Source Term 5 has been Written for JL = co e , w = w e 

a A 

• Eqs. ( 18) - ( 23 ) , (16) and (17) are Solved Using a Central-Difference, Finite-Volume 
Scheme Using Four-Stage Runge Kutta Time Stepping with Added Second- and Fourth- 
Order Dissipation Terms. 

• The Computer Program is a Three-Dimensional Program. 
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• Local-Conical Flow Problem 


• If the Conical Coordinates are Used to Transform the Relative Motion Equation, 
Eqs. (18)-(22), the Resulting Equations will not Represent a Conical Flow. 


® If the Conical Coordinates are Used to Transform the Absolute Motion Equations, 

Eqs. (l)-{7), the Resulting Equations will Represent a Conical Flow for the Steady 
Flow. For the Unsteady Flow, the Problem is Made "Locally Conical" if it is Solved 
at a Fixed Axial Location. 


• In the Three-Dimensional Program, Local Conical Flow Solutions are Obtained at x=l 
by Equating the Absolute Motion in the First and Third Planes: 


(pV) 1tl - (pV), - (pV r ) tll - (pV r ) = (pV r >. + u> x( Pj f t - P(tl ? 1tl ) 


(24) 


<P E) 1t l ■ ((«), * <P E r ) 1± l ' (P E r >, + < 5 r 1 + ? ( - < 25) 


where i = 2. 
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1. Initial Condition: 


For Constant Rolling Problem (w = -w e ),_the Flow Corresponds to a Uniform 
Translation Plus a Rigid Body Rotation -ijxr. 


For Rolling Oscillation (w = -w sin kt e ), the Flow Corresponds to a Uniform 
Translation Only. 


2. Boundary Conditions : 


Normal Momentum Equation is Used on the Wing Surface 
({Sr (v r • n ) = 0, = 0) 


pV . (V • vn) = — + p n • [2wxV + wxr + ojx(o)xr)^ (26) 


In the Farfield, a Uniform Translation Plus the Corresponding Rotation are 
Imposed Outside of the Bow Shock. 
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SYMMETRIC COMICAL FLOW 


This is a verification test case for the three-dimensional program which has 
been solved earlier by using a conical flow program^. Figure 2 shows the results 
for a flat plate sharp-edged delta wing at M m = '2, a = 10° and e (sweep angle) 

using a modified Joukowski transformation^ of 128x64 cells around and normal to the 
wing has been used for the whole computational region. Figure 2 (a--d) shows the 
surface pressure, cross-flow velocity, cross-flow Mach contours and static pressure 
contours. It is clear that two symmetric leading-edge vortices have been captured on 
the suction side along with a weak cross-flow shock under each vortex. The outer bow 
shock is clearly visible in the lower portions of the cross-flow Mach and static 
pressure contours. The surface pressure, Mach contours and static-pressure contours 
match those obtained by using the conical flow program for half of the computational 
region. 



a. surface pressure b. cross-flow velocity 

Fig., 2 Steady symmetric flow around a delta 
wing, M oo =2, a =10°, p=70° 
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THREE-DIMENSIONAL TRANSONIC FLOW 


Figure 3 shows the results of a sharp-edged delta wing for M^ = 0.7, a = 10° and 
aspect ratio of 1.5 using a number of cells of 80x38x48 in the x,°° n and $ direc- 
tions; respectively. The results of Figures 3a and 3b show, from left to right, 
the surface pressure, the static pressure contours and the cross-flow velocity; each 
at the chord stations of 0.52 and 0.81. Figure 3c shows the static pressure 
contours and the cross-flow velocity at the chord station of 1.01 and Figure 3a 
shows the static pressure contours and the cross-flow velocity at the chord station 
of 1.25. At x = 0.81, comparisons of the computed surface pressure with the 

experimental data^ shows that the location of vortex core is well predicted, while 
the value of peak suction pressure under the vortex core is slightly under- 
predicted. At x = 1.25, the static pressure contours and the cross-flow velocity 
show the formation of trailing-edge vortex core and its interaction with the leading- 
edge vortex core. Our cross-flow planes are taken normal to the wing surface. 



a. X =0.52 


b. X =0.01 


Figure 3 . Three-Dimensional Transonic Flow, Standard Euler Set, Sharp-edged delta 
wing, 80X38X48 cell, M^O.7, o»15°, AIM. 5, e 2 =0.12, t^O.QOS, 

(».,b.) 1. Surface Pressure, 2. Static Pressure Contours, 3. Crossflow 

velocity 

Static pressure Contours, 2. Crossflow velocity 
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Figure 3 . Three-Dimensional Transonic Flow, Standard Euler Set, Sharp-edged delta 
wing, 80X38X48 cell, M^O.7, o*15°, AR*1.5, e 2 =0.12. e 4 =0.005. 

(a, b.) 1. Surface Pressure, 2. Static Pressure Contours, 3. Crossflow 

velocity 

(c.,d.) i. static pressure Contours, 2. Crossflow velocity 
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THREE-DIMENSIONAL LOW-SPEED FLOW 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Figure 4 shows the results of a sharp-edged delta wing for M^ = 0.3, a = 20.5° 
and aspect ratio of 1 using a number of cells of 80x38x48 in the x, n> and ? 
directions; respectively. The results of Figures 4a and 4b show, from left to 
right, the surface pressure and the experimental data of Hummel^ the static-pressure- 
coefficient contours and the cross-flow velocity at the chord stations of 0.52 and 
0.81. The results of Figures 4c and,4d show, from left to right, the experimental 

static-pressure-coefficient contours^ (cross-flow planes are normal to wind 
direction), the computed static-pressure-coefficient (cross-flow planes are normal to 
wing surface) and the computed cross-flow velocity at the chord stations 1.02 and 


1.25. 


The results show that the location of leading-edge vortex core is accurately 
predicted, the suction pressure peak is accurately predicted at x = 0.52 but it is 
overpredicted at x = 0.81, the location of the trai 1 ing-edged vortex core is slightly 
off that of the experimental data. A finer grid than the one used in this example is 
expected to give a more accurate prediction. 




Figure 4 . Three-dimensional subsonic flow, isentropic Euler set, sharp-edged 
delta wing, 80x38x48 cell, M. * 0.3, a = 20.5°, AR = 1, e 2 = 

0 . 12 , £4 * 0.005, (a.,b.) 1 . surface pressure, 2 . static presure 
contours, Cp, 3. crossflow velocity 
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UNIFORM ROLLING IN A CONICAL FLOW 

Figure 5 shows the results for a flat plate sharp-edged delta wing which is 
undergoing uniform rolling in the counter-clockwise direction around its axis ox' at 

a constant angular speed o> = = 0.5; where w and to are the dimensional and 

L oo 

dimensionaless angular speeds, i is the wing root chord and U is the freestream 

oo 

speed. The wing angle of attack a = 0 and hence the flow is steady in the rotating 
frame of reference. Figure 5a shows the upper "l" and lower "a" surface 
pressure. As symmetric surface pressure. Figures 5b, 5c and 5d show the 
corresponding cross-flow velocity, cross-flow Mach contours an static pressure 
contours. On the cross-flow Mach contours and under the anti -symmetric leading-edge 
vortices, one notices a weak cross-flow shock. 




b. cross-flow velocity 


Fig. 5 . Uniform rolling of a delta wing, 
M =2, o=0° , (3=70°, uj=0.5 

CD 
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ROLLING OSCILLATION IN A LOCALLY-CONICAL FLOW 


The wing is given a rolling sinusoidal oscillation of u and 6 of w = - w cos kt 

e , and e = - e sin kt, where e mav = u /k and k is the dimensionless reduced 
^ nia a max o 

frequency of oscillation (K = jj— ^ - is the dimensional frequency). Choosing 

0 max = u ^ anc * % = 0*35, the corresponding k = 1.337 and the period of oscillation 

t = 4.699. Figure 6 shows the rolling oscillation motion. 



Point No. Time 9° Sense Figure 


1 

1.07 

-14.85 

CCW 

6 

2 

2.16 

-3.76 

CW 

7 

3 

3.19 

13.5, 

CW 

8 

4 

4.31 

7.46 

CCW 

9 

5 

5.35 

-11.46 

CCW 

10 

6 

6.46 

-10.63 

CW 

11 

7 

7.48 

8.17 

CW 

12 

8 

8.59 

13 

CCW 

13.1 

9 

9.69 

-5.69 

CCW 

13.2 

10 

10.7 

-14.79 

CW 

13.3 

11 

11.8 

1.03 

CW 

13.4 

12 

12.8 

14.8 

CW 

13.5 

13 

13.9 

3.93 

CCW 

13.6 

14 

15.0 

-14.01 

CCW 

13.7 


Fig. 6. Roll angle, angular speed and angular 

acceleration of the rolling oscillation 
motion. 
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ROLLING OSCILLATION (CONTINUED) 


Figure 7 shows the results for the time range t = 0 -1.07. By the end of this 

time, the wing has rolled through an angle e = -14.85° which corresponds to the end 

of the first quarter of the cycle. At t = 0, |oj| = oj 0 and |w| decreases in the 

counter-clockwise direction within that time. Figure 7a shows the surface pressure 
after each 400 time steps covering a total of 2000 time steps. On the upper surface, 
the suction pressure on the left is higher than that on the right, and the suction 
peak is moving in the spanwise direction in the positive z direction. On the lower 
surface, the surface pressure is decreasing on the left side while it is increasing 
on the right side. Figures 7b, 7c and 7d show the cross-flow velocity, cross-flow 
Mach contours and static pressure contours at t = 1.07 and e = -14.85°. At this 
moment, a large leading-edge vortex appears on the left and a small leading-edge 
vortex appears on the right. The cross-flow Mach contours show shocks above and 

below the left leading-edge vortex. It also shows the outer bow shock with varying 

strength (lower portion of the figure). 




a. surface pressure 


b. cross-flow velocity 


Fig. 7 . Rolling oscillation of a delta wing, 
M^, o=10°, P=70° , a>=0. 35, k = 1.337, 

0 max =15 ° • t=0_1 - 07 » 6“0-(-1485 0 ) 
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cross-flow Mach 


d. static pressure 


Fig. 7 . Rolling oscillation of a delta wing, 
M o> =2, a*10°, 0=70° , m- 0.35, k = 1.337, 

6 =15°, t-0-1.07, 0-O-(-14S5°) 




ROLLING OSCILLATION (CONTINUED) 

Figure 8 shows the results for the time range t = 1.07 + - 2.16. Within this 
time, the wing has reversed it direction |u|, and by the end of this time, the 
wing roll angle e = -3.76°. Figure 8a shows the surface pressure after each 400 time 
steps covering the range of time steps from 2,001-4,000. The peak suction pressure 
on the left is decreasing corresponding to a decrease in size of the left vortex 
while that on the right is increasing corresponding to an increase in size of the 
right vortex. Figures8b, 8c and 8d show the cross-flow velocity, cross-flow 
Mach contours and static-pressure contours at t = 2.16, e = -3.76°. 




b. cross-flow velocity 


Fig. 8 . Rolling oscillation of a delta wing, 

M <jo =2 , o=10° , p®70° , (0=0.35, k=1.337, 

6*=15°. t=1.07 + -2.16, e=H4.85°)-(-3.76°) 

ino a 


ORIGINAL PA3£ 

OF POOR QLALr/Y 


280 








ROLLING OSCILLATION (CONTINUED) 

Figure 9 shows the results for the time range t = 2.16 + - 3.19 during which the 
wing is rotating in the clockwise direction with decreasing |u|. By the end of this 
time, the roll angle is e = 13.5°. Figure 9a shows the surface pressure covering 
the range of time steps 4,001-6,000. The peak suction pressure on the left is moving 
to the left as the vortex is disappearing, and an attached flow is forming. Thepeak 
suction pressure on the right is moving inboards to the left, while the shock under 
the vortex is growing. Figures 9b, 9c and 9d show the cross-flow velocity, cross- 
flow Mach contours and static pressure contours at t = 3.19 and e = 13.5°. 


0 

-0 r«.?T7M 



a. surface pressure b. cross-flow velocity 

Fig. 9. Rolling oscillation of a delta wing, 

M oq = 2, a=10°, P=70° , to=0.35, k-1.337, 

0* =15°, t=2.16 + -3.19, 0 = (-3.76° )- 

(+13.6°) 
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c. cross-flow Mach 


d. static pressure 


Fig. 9.. Rolling oscillation of 
M =2, a=10° , 0=70°, u)=0 

(+13.5°) 


a delta wing, 
.35, k=1.337, 

0=(-3.76° )- 





ROLLING OSCILLATION (CONTINUED) 

Figure. 10(a-d) shows the results for the time range t = 3.19 + - 4.31 during 
which the winq has reversed the direction of rotation from the CW to CCW, and |to| has 
reached a zero value and then increases. The vortex on the left is growing, and the 
corresponding peak suction pressure is increasing and moving inboards to the right. 
The vortex on the right is flattening, and its peak suction pressure is decreasing 
and movino inboards to the left. By the end this time t = 4.31; the roll angle 
9 = 7.46°I 


-.4 

0 !=.343E*1 

0 rc.3B7M 
O !=-388E*l -* 3 

r °4v^ 

ilk l 

i i i i P 

_ i i i i ; 

0 -.8 -.6 -.4 -.2 

.4 

i -2 .4 -6 .0 10 

0 = 746 ° 


Z/B 



a. surface pressure 


b. cross-flow velocity 


Fig. 10. Rolling oscillation of a delta wing, 

14^=2, a=10° , p=70°', w=0.35, k = 1.337, 

0 max =15 °* t= 3-19 + -4.31 , 9 = (+13.5 0 )- (+7.46°) 
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ROLLING OSCILLATION (CONTINUED) 


Within the time range t = 4.31 + - 5.35, the wing is still rotating in the CCW, 
and LI has reached its maximum value and then decreases. At t = 4.7, the wing has 
already completed one cycle of oscillation. The peak suction pressure on the left is 
moving inboards to the right, and the peak suction pressure on the right is moving 
outboards to the right, Figure 11 (a--d). By the end of this time t = 5.35, the roll 
angle e = -11.46°. 
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o r* 


-.5 

-.1 

-.3 







O 

O !t.473£.1 

A It.SMM 
k Tt .53tf ♦! 


,^00 

K. 

A.- 







J L 


i — i 

l 


.8 


Z/B 


if ttf HHfflBM 


eT- 


* l 


0 =-11.46 



a. surface pressure 


b. cross-flow velocity 


Fig. 11. Rolling oscillation of a delta wing, 

14^=2, a=10° , P=70°, w=0.35, k = l . 337 , 

9°° =15°, t-4.31 + -5.35, 0 = (+7.46° )- (-11.46°) 

max 
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c. cross-flow Mach 


d. static pressure 


Fig. 11, Rolling oscillation of a delta wing, 
1^=2, a=10° , (3=70° , u=0.35, k=1.337, 

0* «15°, t=4.31 + -5.35, 0=(+7.45° )“ 

ma X 
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ROLLING OSCILLATION (CONTINUED) 

Figures 1 2a- 1 2d and 13a-13d show the results covering the ranges of time steps 
10,001-12,000 and 12 ,001-14 ,000 , respectively. Figures 5-11 cover a total of 1.56 
cycles of oscillation. The results show the successive increase, decreased and 
motion of the left and right vortices and their corresponding peak suction pressure. 
They also show the formation and disintegration of the shocks below and above the 
vortices, and the motion and strength variation of the outer bow shock. 




b. cross-flow velocity 




Fig. 12. Rolling oscillation of a delta wing, 

M oo =2, a=10° , p=70° , co=0.35 p k = l. 337, 
e” =15°, t=5.35 + -6.46, 0 = (-ll. 46°)- (-10.63°) 

ITla X 
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c. cross-flow Mach d- static pressure 


Fig. 13. Rolling oscillation of a delta wing, 

Vlj=2, a=10° , p=70°, u)=0.35, k=1.337, 

0°° =15°, t=6.46 + -7.48, 0 = (1O.63° )-(+8.17°) 
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ROLLING OSCILLATION (CONCLUDED) 

Steady state oscillation response is reached after 3 cycles of transient 
response, Figure 14 (a— g). 



Fig. 14 . Rolling oscillations from t = 7.48 
to t = 15.0, steady state oscilla- 
tion is reached. 
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Fig. 14. Continued. 
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Continued. 
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ROLLING OSCILLATION, LIFT AND ROLLING-MOMENT COEFFICIENTS 

Figure 15 shows the time history of the lift and rolling-moment coefficients 
along with roll angle variation. Steady-state oscillation response is reached by 
the third cycle, while the phase angle between C^ and ©is 90°, C M is in phase with o>. 
Although C M and e have the same frequency, C^ snows twice the value of that frequency. 

.10 
.oe 
.06 
.04 
.02 

°U 
■02 
.04 
.06 
.06 

■ 3b 1 I I I I I I l-.JO 

0 I 2 3 4 S 6 7 

T*K/PI 

FIGURE 15. TIME HISTORY OF THE LIFT AND ROLLING-MOMENT COEFFICIENTS 
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3. Integral Solution of Full -Potential Equation 
With and Without Embedded Euler Domains 


Formulation 


Full Potential Equation (Shock Capturing, SC; Shock Capturing-Shock Fitting, SCSF) 


$ + <$ =» G 

XX yy 


(1) 

G " nr <p x \ * 

P y v 

(2) 

P ■ (1 + M 2 

2 <* 

(1 - » 2 - t 2 )] 1 ^ 1 
x y J 
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V<I> • n ** 0 

on g(x, y) =» 0 

(4) 

V<J> e 

00 

away from g 

(5) 

A C 1 - 0 

p l TE 


(6) 


Integral Solution of Velocity Field With Explicit Shock Surface Contribution 
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Shock Fitting 
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Pressure Coefficient 
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IW 


Euler Equations (Integral Equation With Embedded Euler Domains) 


M a? 

St bx by e ® 

(15) 

Q • Ip, pu, pv, pej 6 

(16) 

— 2 

E = ^ u » P u + P» puv, puh]* 

(17) 

«• 2 
F = [ P V » P UV » pv + p, pvh]* 

(18) 

( Y +l)p + (u + V )/2, h = e + p/p 

(19) 
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Method of Solution 


• Shock-Capturing Shock-Fitting (SCF) Scheme: 

Shock Capturing Part 

1. In Eq. (7), Set G = q s = 0, Use Eq. (7) to Satisfy B.C.S; Eqs. (4), (6). 
Find qg and y g (piecewise linear distribution). 

2 . Calculate Initial Values of G at the Centroids of Field Elements 

(G = M* u x ). 

3. With q $ = 0, Eq. (7) is Used to Satisfy B.C.S. 

4. Calculate p and G Using Eas. (2), (3) (Type Finite Difference is Used 

for p x * p y) * 

5. Steps (3) and (4) are Repeated Until Captured Shock Location is Fixed. 
Shock Fitting Part 

6. Shock Panels are Introduced, q $ - term in Eq. (7) is activated. 

7. Eqs. (8) and (12) are Used to Find q s and p, and Eqs. (9)-(12) are Used to 

Cross the Shock. 

8. Iterative Procedure is Continued Until Convergence. 
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• Integral Equation With Embedded Euler (IEEE) Scheme (For Strong Shocks) 

1. Shock Capturing is Used to Locate the Shock. 

2. A Computational Fine-Grid Domain is Constructed Around the Shock for Euler 
Computations. 

3. With B.C. and I.C. Found From Step (1), Euler Eqs. ( 15)-( 19) are solved by Using 
a Central-Difference Finite-Volume Solver with Four-Stage Runge-Kutta Time 
Stepping and Added Second- and Fourth-Order Dissipation. 

4. Fixing q Values Found From the Euler Calculations, the Integral Ea. is Used to 
Update the B.C. 

5. The Iterative Procedure is Continued Until Convergence is Achieved. 



Fig. 16 Integral Equation Grid with an 
Embedded-Euler Domain. 
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SHOCK-FREE FLOW 


The first step to validate the computer program is to check the sensitivity of 
the IE solution to the size of the computational domain. Figure 17 shows the 
solutions for the NACA 0012 airfoil at = 0.72 and a = 0° using vortex panels only 
on the airfoil surface. We used a total°°of 140 vortex panels on the airfoil surface 
and a 64x60 field elements around the airfoil. The solutions show the surface 
pressure using two sizes of the computational domains; 2x1.5 and 3x2.5. In Figure 18, 
we repeat the same test for a lifting case of the same airfoil at M = 0.63 and 
a = 2°. The results of these two cases show that a computational domain of 2x1.5 
gives as accurate solutions as those of the 3x2.5 computational domain. 



Fig. 17 Effect of the Computational Domain Fig. 18 Effect of the Computational Domain 
Size, Surface Vortex Panels, NACA Size, surface Vortex Panels, NACA 

0012, M a =0.72, <x=0°. 0012, M =0.63, a=2° . 
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SHOCK-FREE FLOW 


The second numerical test is aimed at comparing the results of the IE solution 
using vortex panels only and source panels only with the solution of Euler 
equations^. Figure 19 shows the results of this test for the NACA 0012 airfoil at 
= 0.7 and a = 0°. The computational domain is 2x1.5, and the same numbers of 
surface panels and field elements as those of Figure 17 have been used. It is clear 
that the IE solution with surface vortex panels is superior to that of the source 
panels. 



o .2 .4 .6 .8 i.o 
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Fig. 19 Comparisons of IE Solution with 

Surface Vortex Panels and Surface 

Source Panels with Euler Solution, 

NACA 0012, M =0.72, ot=0° . 

00 
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TRANSONIC FLOW 


First, we present a numerical test case to show the effect of introducing the 
shock panels and their fitting as explained earlier. Figure 20 shows a comparison 
between the shock capturing results and the SCSF-scheme results for the NACA 0012 
airfoil at = 0.8 and a = 0°. It is clear that the SCSF-scheme sharpens the 

shock, as expected, with this relatively coarse grid. Next, we compare the SCSF- 
scheme with the experimental data and other computational results. Figure 21 shows 
the results of the SCSF-scheme for NACA 0012, M = 0.8 and a = 0°, along with 

comparisons with the computational results of Garabedian, Korn and Jameson 0 and the 
experimental data taken from reference 7. The SCSF-scheme took 12 iteration cycles 
of shock capturing (SC) and 13 cycles of shock fitting (SF) to achieve convergence. 
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TRANSONIC FLOW 


Figure 22 shows the results of the IFFF-scheme for the same case along with a 

comparison with the computational results of Jameson®, who also used the finite- 
volume Fuler scheme with four-stage Runge-Kutta time stepping. In the present I EEE- 
schene, the embedded Euler domain has a size of 0.5x0. 6 around the shock region with 
a grid of 25x30. This case took 10 iteration cycles of SC, 250 times cycles of Euler 
iterations to achieve a residual error of 10"'* and 5 IE cycles to update the Euler 

domain boundary conditions. 



0.0 .2 .4 .6 .8 1.0 
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Fig. 22 Integral Equation with Embedded- 
Euler Domain Solution, NACA 0012, 
M =0.8, a=0° . 
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TRANSONIC FLOW 


Figures 23 and 24 show the results of the SCSF-and IEEE-schemes for NACA 64A010A, 
M = 0.796, a = 0° along with comparisons with the computational results of Edwards, 
Bland and Seidel 9 who used the TSP-equation, and the experimental data taken from 
reference 9. With the SCSF-scheme, the numbers of SC and SF iteration cycles to 
achieve convergence are the same as those of the case presented in Figure' 21. With 
the IEEE-scheme, the embedded Euler domain has a size of 0.7x0. 6 with a grid size of 
35x30. This case. Fig. 24, took 10 iteration cycles of SC, 130 time cycles of Euler 

iterations to achieve a residual error of 10”^ and 3 IE cycles to update the Euler 
domain boundary conditions. 



Fig. 23 Integral Equation Solution with Fig. 24 Integral Equation with Embedded- 

SCSF-Scheme, NACA 64A010A, Euler Domain Solution, NACA 64A010A, 

M =0.796, a=0° . M^O.796, a=0° . 
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TRANSONIC FLOW 


Figures 25 and 26 show the results of the SCSF- and IEEE-schemes for the lifting 
case of NACA 0012, = 0.75 and a = 2° along with the computational results of 

Steger and Lomax*®, and the experimental data taken from the same reference. The 
size of the grids and the number of iteration cycles used to achieve convergence are 
the same as those of the cases given in Figures 21 and 22. 




Fig. 25 Integral Equation Solution with 
SCSF-Scheme, NACA 0012, M =0.75, 
a=2° . 


Fig. 26 Integral Equation with Embedded- 
Euler Domain Solution, NACA 0012, 
M^O.75, a =2° . 
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TRANSONIC FLOW WITH STRONG SHOCKS 


For stronger shocks than those considered above the IE computational domain is 
extended in the longitudinal and lateral directions and so is the embedded Euler 
computational domain. The Euler domain is extended beyond the trailing edge to allow 
for the vorticity to be shed downstream where the overlapping region with the IE 
equation exists. Figure 27 shows a typical computational domain with details of the 
embedded Fuler domain. 



Fig. 27 Embedded Euler Domain and Grid for 
Strong Shocks. 


! 
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TRANSONIC FLOW WITH STRONG SHOCKS 


Fioure 28 shows the results of the IEEE for NACA 0012, M = 0.812 and a = 0° 

along with the experimental data of reference 7. In Figure 29, the results of the 
IEEE for NACA 0012, fl = 0.82 and a = 0° are shown along with the three-dimensional 
solution at the wing root chord of Tseng and florino , who use the IE for the TSP, 
and the results of reference 11. The size of the embedded Euler domain for these 
cases is 0.8x0. 8 and its grid size is 40x40. 
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Fig. .28 Integral Equation with Embedded- 
Euler Domain Solution, NACA 0012, 
M =0.812, a=0° . 
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TRANSOM I C FLOW WITH STRONG SHOCKS 


Figure 30 shows the results of the IEEE for NACA 0012, M m = 0.84 and a = 0° 
along with comparisons with the nonisentropic FP-solution of Whitlow, et al and 
the Euler equations solution of Jameson®. The size of the embedded Euler domain for 
this case is 1. 5x1.0 and its grid size is 60x40. This case took 10 IE iteration, 300 
time cycles of Euler iterations and 3 IE cycles to update the Euler domain boundary 
conditions. 



Fig. 30 Integral Equation with Embedded- 
Euler Domain Solution, NACA 0012, 
M =0.84, <x=0°. 

CD 
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4. Concluding Remarks 


• Two Methods have been Presented for Inviscid Transonic Flows: 

- Unsteady Euler Equations in a Rotating Frame of Reference for Transonic-Vortex 
FI ows. 

- Integral Solution of Full -Potential Equation with and without Embedded Euler 
Domains for Transonic Airfoil Flows. 

• The Computational Results Covered: 

- Steady and Unsteady Conical Vortex Flows 

- Three-Dimensional Steady Transonic Vortex Flow 

- Transonic Airfoil Flows 

• The Results are in good agreement with Other Computational Results and 
Experimental Data. 

• The Rotating Frame of Reference Solution is Potentially Efficient as Compared 
with the Space-Fixed Reference Formulation with Dynamic Gridding. 

• The Integral Equation Solution with Embedded Euler Domain is Computationally 
Efficient and as Accurate as the Euler Equations. 

• Currently the Rotating Frame of Reference Eule Solver is Applied to Three- 
Dimensional Unsteady Transonic-Vortex Flows. The IEEE-Scheme is Being Extended 
to the Unsteady Transonic Airfoil Calculations. 
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